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ABSTRACT 

Stellar activity signatures such as spots and plage can significantly limit the search for extrasolar planets. Current mod- 
els of activity- induced radial velocity (RV) signals focused on the impact of temperature contrast in spots predicting the 
signal to diminish toward longer wavelengths. Contrary to this is the Zeeman effect on radial velocity measurements: 
the relative importance of the Zeeman effect on RV measurements should grow with wavelength because the Zeeman 
displacement itself grows with A, and because a magnetic and cool spot contributes more to the total flux at longer 
wavelengths. In this paper, we model the impact of active regions on stellar RV measurements including both tem- 
perature contrast in spots and line broadening by the Zeeman effect. We calculate stellar line profiles using polarized 
radiative transfer models including atomic and molecular Zeeman splitting over large wavelength regions from 0.5 to 
2.3 nm. Our results show that the amplitude of the RV signal caused by the Zeeman effect alone can be comparable 
to that caused by temperature contrast; a spot magnetic field of ^1000 G can produce a similar RV amplitude as 
a spot temperature contrast of ^lOOOK. Furthermore, the RV signal caused by cool and magnetic spots increases 
with wavelength contrary to the expectation from temperature contrast alone. We also calculate the RV signal due to 
variations in average magnetic field strength from one observation to the next, for example due to a magnetic cycle, 
but find it unlikely that this can significantly influence the search for extrasolar planets. As an example, we derive the 
RV amplitude of the active M dwarf AD Leo as a function of wavelength using data from the HARPS spectrograph. 
Across this limited wavelength range, the RV signal does not diminish at longer wavelengths but shows evidence for the 
opposite behavior consistent with a strong influence of the Zeeman effect. We conclude that the RV signal of active stars 
does not vanish at longer wavelength but sensitively depends on the combination of spot temperature and magnetic 
held; in active low-mass stars, it is even likely to grow with wavelength. 
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1. Introduction in the shape of stellar line profiles on timescales similar to 

planetary orbits can be ca used by several mechanisms, e.g. , 
The precise determination of radial velocities (RV) and by the transit of a planet (|Rossiterlll924l: lMcLaughIinlll924[) 
their temporal variations is a key data analysis method in or stellar activity. The latter poses a number of problems 
stellar astrophysics. It is applied to detect extra-solar plan- to stellar radial velocity measurements: first, active stars 
^ ets and to measure their projected mass, which requires are typically fast rotators implying wider line profiles re- 
^ ' high p recision RV data pushing to the order o f ms"-^ and ducing achieva ble accuracy in a gi ven radial velocity mea- 
i below (jMavor fc U drvl 120081: lUdrv et al.ll2009D . Radial ve- surement (;e.g.. iBouchv et al.|[200l[) . Second, activity is be- 
. !^ ' locities can also be used to determine periodic motion of lieved to be caused by magnetic areas that can produce 
K> ] the stellar surface enabling asteroseismology to reveal sen- cool spots or hot plage and are in general not symnietri- 
^ , sitive information on fundamental stel lar parameters in- cally distributed over the star; this allows the reconstruc- 
ts ■ eluding a vie w into the interior of stars (jCunha et al.ll2007t tion of surface m aps in tomographic imaging studies (e.g., 
iBazot et al.l l20TTI). iVoet fc Penr^ 11983). Third, magnetic regions can sup- 
Measurement of RV time series allows the detection of P^ss convective motion and alter the signature of stellar 
the projected motion of a star or its surface from their spec- convective blueshifts (see Gray 2009; Mcunicr et al. 201Q). 
tral lines. For the detection of planets, the assumption is ^ool spots co-rotating with the stellar surface are well stud- 
that the shape of a spectral lines does not vary with time for the Sun and much larger spots are known to exist 
so that its centroid position provides information about the on other stars (|Berdyugina| [20051; [Strassmeieij 120091). Line 
projected vefocity of the star. The relative radial vefocity Profile distortions due to these features can lead to signif- 
shift between two epochs is measured either by searching for icant shifts of the line barycenter introducing an offset to 
the best agreement between two spectra with radial velocity ^^eal central line position (e.g.,|Saar et al.||1998D. Useful 
as a free parameter, or by locating the centroid position in information for the correction of stellar RV curves for activ- 
a cross correlation profile calculated from the spectrum and signatures can be provided by activity indicators hke the 
some template. Both methods assume that the spectral line strength of chromosphcric emission lines or absorption line 
shape is identical. It is well known, however, that variations bisectors ( [Santos et al.,,2000 ,). Several examples exist where 
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periodicities in radial velocity curves were interpreted as 
due to planetary companions but that later were retracted 
since the reason for RV variability was found to be stel- 
lar activity. Additionally, differences in RV amplitude be- 
tween optical and infrared bands have been detected in a 
few systems (e g.. iMarti'n et al.l l2006l: iHuelamo et"alll2008l: 



iMahmud et al.ll201l[ ). 

The effect of cool active regions on radial velocity 
measurements due to temperature co ntrast alone has 
been investigated in detail by, e.g .. ISaar fc Donahue! 
(I1997D: Dcsort et al. 



detail by, 

(2007); ' Ma karov et al.l 



(I2009D. 



iLagrange et al.l ((20100; IReiners et al.l ( 20101 ): Barnes et ahl 
(|2nilD . Magnetic flux tubes on the Sun are known to 
exist also in the so-ca lled network and plage areas (e.g., 
ISolanki fc Stenfldll984D . These regions show high magnetic 
fields well above B = 1 kG, temperature contrasts of a 
few hundred K, and occupy a much larger fraction of the 
solar surface than cool spots do. Due to their relatively 
low temperature contrast but large filling factor, plage are 
difficult to characterize on surfaces of stars other than the 
Sun. The impact of bright regions on RV measurements 
(together with the impact of inhibi ted convection) was 
investigated by iMeunier et al.l (|2010( ) using solar 2D im- 
ages as template for the spatial distribution of bright and 
dark areas. In principle, the effect of bright areas on RV 
measurements is comparable to the one of cool spots, but 
since flux contrast is reversed, the RV signal has opposite 
sign. RV distortions from adjacent bright and dark areas 
can therefore partially cancel out. The total contribution 
of plage to the variability of solar irradiance is larg er tha. n 
the contribution of cool spots (e.g.. Frohlich &: LeanillQQSl) . 
but the two so far cannot be distinguished in any other 
star. Therefore, plage and quiet stellar regions are often 
described by one single atmospheric component, which 
probably is somewhat hotter than the "quiet" atmosphere. 
Line profile distortions in stars other than the Sun are 
often described in a two-temperature model defining a cool 
spot component in active stars. The real effect is likely a 
result of spot and plage variability. 

In general, co-rotating active regions can lead to spu- 
rious radial velocity variations in phase with the stellar 
rotation period. The amplitude of the variation depends 
on the projected rotational velocity of the star, v sini, the 
fractional surface coverage of the spotted area, /, the tem- 
perature contrast between "quiet" star and active regions, 
and the wavelength used for observations. Since the (abso- 
lute) flux ratio between active regions and photosphere is 
smaller at longer wavelengths (independent of whether the 
active region is cool or hot), it is expected that radial veloc- 
ity distortions due to stellar activity are lower at infrared 
than at optical wavelengths. It has therefore been claimed 
that radial velocity-based planet searches in very cool stars 
and brown dwarfs (these objects are typically very ac- 
tive) should be optimally co nducted at infrared wavelengths 
rather than at optical fe.g.. lMartm et a'Lll2006D . 

The aim of this paper is to investigate the impact of 
magnetic fields on the spectral appearance of active stars. 
It is not only the temperature difference that affects stellar 
line profiles, but also the field itself that may introduce sig- 
nificant distortions through the Zeeman effect. The Zeeman 
effect operates opposite to temperature contrast with more 
significant influence at longer wavelengths. Our aim is to 
model the radial velocity signal of active regions including 
both temperature and magnetic effects. After introducing 



our model techniques and demonstrating the general pic- 
ture in a few toy model cases, we conclude that the Zeeman 
effects likely plays a very signiflcant role in the determina- 
tion of radial velocities in active stars. 



2. Zeeman splitting in stellar spectra 

Our picture of stellar surface properties is motivated by the 
appearance of the Sun, where rising magnetic flux tubes in- 
hibit convection in active regions and produce both hot and 
cool surface regions. Strongest flux concentrations are ob- 
served in cool spots with typical temperatures several 100- 
1000 K below photospheric average and typical magnetic 
flux densities of several 100-1000 G (Solanki 2003,). 

The general framework of stellar active regions is con- 
sistent with observations from more active stars and stars 
of different spectral type ()Strassmeieij|2009h . An important 
ingredient is that active regions differ from the quiet pho- 
tosphere in both temperature and magnetic properties; in 
particular, large starspots are believed to be both cool and 
magnetic. Owing to our limited ability to measure local- 
ized magnetic flelds in other stars, however, we lack good 
understanding of magnetic fields in stars other than the 
Sun. In particular, we have no empirical evidence of the 
relation between magnetic field and spot temperature in 
very active stars. Nevertheless, we may find it reasonable 
to assume that local magnetic field strengths in other stars 
are on the same order as in the Sun, and that large active 
regions can be similar to very large sunspots. Empirical 
results on average magnetic fields in very active low-mass 
stars are consistent with this picture finding very strong 
average surface fiel ds at the kilo-Gauss le vel in mid- and 
late-M dwarfs (e.g.. iReiners fc Basril[2007t ). 

The appearance of a spectral line in the presence of a 
magnetic field is determined by the Zeeman effect: each en- 
ergy level with a total angular momentum quantum num- 
ber J splits into (2 J -I- 1) states of energy with differ- 
ent magnetic quantum numbers M. In absence of a mag- 
netic field, the transition energy is unique but it splits into 
three groups of transitions according to the change in the 
magnetic quantum number M invoked by the transition 
(AM = —1,0,-1-1). The appearance of the spectral line also 
depends on the geometry of the field, but this effect is often 
neglected assuming an "homogeneous" distribution of fi eld 
lines over the stellar surface. We refer to iReiner 1 (I2012D for 
a deeper discussion of magnetic field observations. 

For the context of radial velocity measurements, we are 
interested in the amplitude of spectral line deformations 
caused by magnetic fields. The two cr-groups with magnetic 
quantum numbers M = —1, +1 are shifted with respect to 
the TT-group {M — 0) by an amount that depends on the 
level's quantum numbers, condensed in the Lande-factor 
g, and is proportional to the magnetic field B. In velocity 
units, the displacement can be equated as 



lAgXB, 



(1) 



with V in ms~^, A in /um, and B in Gauss. The Lande-factor 
g is of order unity. Eq.[T] has two important implications: 
1) The typical displacement of Zeeman components in the 
presence of magnetic fields is on the order of lms~^ G~^. 
For typical field strengths of solar active regions (100- 
1000 G), fiux from active regions can therefore be displaced, 
through the Zeeman effect, by several hundred ms~^; 2) 
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in velocity units, the displacement is proportional to wave- 
length, A, of the spectral line under consideration. Thus, the 
displacement is larger for longer wavelengths, which is con- 
trar y to the displacemen t from temperature contrast (see, 
e.g.. iReiners et alll201Cl[ ). We note that RV signatures due 
to temperature contrast always diminish with wavelength 
independent of whether they are hotter or cooler than the 
rest of the star. The translation from the displacement of 
Zeeman components, Awzccman, into a shift of the line pro- 
file barycenter, Av^rad, is non-trivial and subject of this 
paper. The typical amplitude of Awzocman together with 
typical field strengths on stars show that even if the net sig- 
nal (Awvrad) in & spectral line only would be a few percent 
of Auzeeman, this Still would casily be in the range relevant 
for detecting planetary orbits through the RV method. 

For the purpose of this paper, we are interested in the 
effect of Zeeman splitting on stellar radial velocity curves. 
A constant average magnetic field may affect the overall 
shape of a line profile with respect to the non-magnetic 
case, but as long as this profile is not time- variable, it is 
not relevant for radial velocity analysis. On the other hand, 
any time-variability in the magnetic properties of stellar 
surfaces can have significant consequences for stellar ra- 
dial velocity curves. Variability can be caused either by 
localized magnetism on time-variable (projected) surface 
areas, like co-rotating spots, or by intrinsic time- variability 
of the magnetic field observed on the visible hemisphere 
(e.g., magnetic cycles). We will consider both types of line 
profile variability in the following. 

3. Co-rotating magnetic spots 

Spectral lines of rotating stars are broader than lines of non- 
rotating stars because of the Doppler effect: light emerging 
from the area of the star rotating towards us is blue-shifted 
while light coming from the area moving away from us is 
red-shifted. The net effect is a characteristic spectral line 
broadening allowing precise measurement of the star's pro- 
jected rotational velocity, v sinz (e.g., iGray 2005). If a cer- 
tain area of a star is cooler, this area emits less flux than 
other regions. In absorption lines, this leads to a character- 
istic bump at the position in the line profile corresponding 
to the local velocity projected towards the observer. In a 
similar fashion, the flux emitted from co-rotating magnetic 
spots will alter the appearance of a Doppler-broadened line 
profile, and the shape of the overall profile will change when 
the star is observed at different rotational phases. 

Stellar convective blueshift is another result of Doppler 
broadening in spectral lines. Due to the imbalance between 
rising, hot plasma and cooler downflows, s pectral lines are 
generally blue-shifted in sun- like stars fe.g..'Gravl l200l . In 
magnetic areas, convective motion is suppresse d, wh i ch can 
result in significant RV signatures. iMeunier et al.l (|2010[ ) 
have modeled this effect assuming that in active regions 
the average convective blueshift seen in all lines is attenu- 
ated by Av = 190 ms~^ perpendicular to the solar surface. 
IMeunier et a l. (2010) find that in the Sun, the RV signature 
from convective blueshifts dominates the activity-induced 
stellar RV signal; with an amplitude of several ms"""^, the 
convective RV signal is larger than the signal due to fiux 
contrast by a factor of a few. In our study, we do not take 
into account the signal from convective blueshifts because 
we are concentratin g on the addit i onal effect of Zeeman 
splitting. Following Meunier et al.l (|2010l) . one can argue 
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Fig. 1. Sketch of the effect of a magnetic spot on radial 
velocity measurements. Top: rotating star with large spot, 
covering 5% of the visible surface if observed at disk cen- 
ter, seen at three phases. The star is rotating with the spot 
approaching in the left column, centered in the middle col- 
umn, and rotating out of view in the right coltunn. Second 
row: stellar line profile broadened by rotation but without 
spot flux, showing the impact of a dark star spot. Third 
row: flux emerging from spot only, assuming same temper- 
ature as photosphere. Row four: sum of surrounding plus 
the spot region. Row five: Radial velocity determined from 
fitting a Gaussian profile to the spectral lines in row four. 
In rows three to five, black solid lines show the case for no 
magnetic field (and zero spot temperature contrast). Blue 
and red dashed lines show cases in which the spot area has 
B = 1 kG and 2kG, respectively (A = 1.2 /zm and g = 1.0). 



that the RV variations due to blueshift can be separated 
from the effects of the line shape (fiux contrast and Zeeman 
splitting). 



3.1. Toy model 

In the following, we calculate a line profile for an artificial 
star rotating at v sini = 2 kms""'^ and simulate a line profile 
distortion due to an artificial spot with the same tempera- 
ture as the quiet photosphere. The distortion is only due to 
Zeeman splitting of the line emerging from the spot. In our 
examples, we calculate the rotational phases corresponding 
to maximum displacement of the line center, i.e., we com- 
pute the semi-amplitude K of the app arent radial velocit y 
curve due a magnetic spot (see Fig. 8 of lReiners et al.|[2"010[ ). 
The radial velocity is then calculated from cross-correlating 
the undistorted template spectrum with the spectrum of 
the spotted star. The barycenter of the cross-correlation 
function is determined by fitting it with a Gaussian profile. 
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The top panel in Fig.[T] shows a sketch of rotating star 
at three rotational phases. The stellar surface shows a large 
magnetic spot that is rotating into view in the left column, 
is centered in the middle column, and is rotating out of view 
in the right column. The spot covers 5% of the visible sur- 
face at the center of the disk but it appears smaller if viewed 
closer to the limb. The second (from top) row of Fig.[T]shows 
a stellar line profile broadened by rotation with projected 
rotation velocity v sini = 2kms~^ but without any flux 
coming from the spot region; the total flux is reduced by 5% 
in this (central) example. This shows the typical line profile 
when looking at the temperature effect only. The third row 
shows the flux emerging from the spotted area only, assum- 
ing it has the same temperature as the surrounding photo- 
sphere. In general, large sun- or starspots are believed to 
be cooler than the surrounding atmosphere, but we look at 
the case of identical temperatures first to isolate the impact 
of the Zeeman effect from the temperature contrast. Row 
four shows the sum of the surrounding region (star with- 
out spot) plus the spot region. In panels three and four, 
black solid lines show the case for no magnetic field. Here, 
the profile is just the undistorted spectral line. Dashed blue 
and red lines show cases in which the spot area harbors a 
magnetic field with an average field strength of 1 kG and 
2 kG, respectively. A spectral line at A = 1.2 /im and g — 1.0 
is assumed. The spectral line emerging from the spot region 
is broader and shallower in the magnetic case. In row four, 
the effect of one co-rotating magnetic spot (again, without 
any temperature difference) on a spectral line is displayed. 
Clearly, the Zeeman effect significantly distorts the line pro- 
file and consequently shifts the apparent center of the line. 
Finally, in the bottom panel, we quantify this by showing 
the center of the spectral line as derived from a fit assum- 
ing a Gauss function. We note that the radial velocity shift 
is a consequence of the non-axisymmetric field distribution 
implying that polar spots or other axisymmetric configura- 
tions cannot introduce RV shifts through this mechanism 
(but see SectionlH). We can co mpare Fig.[T] of this work to 
Fig. 8 in iReiners et al.l ()2010l) to see that the line profile 
deformation induced by the Zeeman effect is similar to a 
deformation induced by cool spots. In particular, the radial 
velocity signal from a magnetic spot through the Zeeman 
effect has the same sign as the signal from a cool spot due 
to flux contrast. 

The amplitude of the apparent radial velocity shift in 
our example is approximately 100 ms^^ for B — 2kG, and 
approximately 25ms~^ for B = 1 kG. It is worth noting 
that the projected rotational velocity chosen for this exam- 
ple, V sini — 2kms^^ is similar to solar rotation, which is 
considered rather slow compared to typical late-type stars. 
Nevertheless, the amplitude is significant for the precision 
required for the RV accuracy level needed for planet search. 
For more r apidly rot ating stars the distortion scales with 
V sin I (cD. iReiners et a l. 20101. Because of the similarities 
between temperature and Zeman RV signature s, line profile 
diagnostics like bisectors (jPesort et al.ll2007ll can be use- 
ful tools to investigate activity-related reasons for velocity 
shifts in observed stellar spectra. 

3.2. Dependence on wavelength and field strength 

In Fig. [2] we show the RV amplitude due to the Zeeman 
effect calculated from one spectral line with 5 = 1.0 as a 
function of wavelength. The line distortion is induced by a 
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Fig. 2. Toy model radial velocity signal due to the Zeeman 
effect from a co-rotating magnetic spot with zero temper- 
ature contrast. The signal is calculated for three different 
field strengths inside the spot (green: 100 G, black: 600 G, 
and red: 4000 G); for each case three different spot sizes 
are used (/ = 1 %, 3 %, and 10 %). Grey dashed lines show 
an analytical approximation of the scaling using Eq.[2] The 
dotted line visualizes the 1 ms~^ limit. 

magnetic spot with one out of three different values of mag- 
netic field strength inside the spot {B — 100 G, 600 G, and 
4000 G) and one out of three different spot filling factors 
(/ = 1 %, 3 %, and 10 %). While not much is known about 
the geometric concentration of small magnetic areas on cool 
stars, the total magnetic energy assumed in our examples is 
easily justified by observations of cool star magnetic fields 
einersi |2012[ ) ; our example stars have average fields of 
Bf = 1-400 G (concentrated in one single spot) well in the 
range of average fields observed that can be as strong as 
several kG (Reiners et al. .2009 ). 

In all cases, the apparent radial velocity shift scales with 
spot size and higher field strength introduces larger radial 
velocity shift. For relatively weak fields {B ^ IkG), the 
signal also grows with wavelength as can be expected since 
the Zeeman effect does. The amplitude of the RV signal 
during stellar rotation due to the signature of the spot can 
be approximated by the following scaling relation: 

AWrad.toy = COnSt X / (BA)^, (2) 

with / the relative fraction of the spot area, B the mag- 
netic field inside the spot, and A the wavelength. The grey 
dashed lines in Fig. [5] illustrate the scaling of Eq.jU the ra- 
dial velocity shift is proportional to the filling factor, and 
the dependence on B is identical to the one on A, which is 
consistent with Eq.[TJ 

The amplitude of the radial velocity shift is below 
Ims^^ as long as the field inside this one spot is on the 
order of i? = 100 G. In this case, only relatively high spot 
coverage (10%) at long wavelengths (A > 2 /im) can cause 
signals larger than ■^l ms~^. For the Sun we know that field 
strengths inside a spot can easily be larger than 100 G (but 
note that spots with large fields are typically cool, which 
is not considered in this simple model). In our simulation, 
spots with several hundred Gauss field strengths and filling 
factors of a few percent can introduce radial velocity signals 
well above the ms~^ level. For example, the signal of a spot 
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with / = 3 % and B = 600 G causes a radial velocity signal 
with an amplitude of Awrad ~ Ims"^ at A = 500 nm. The 
signal grows with wavelength up to Auiad ~ 10 ms^^ at 
A = 2200 nm. 

The case of very strong magnetic fields (4000 G) inside 
the spot region shows a somewhat different behavior. While 
the radial velocity shift in this case is larger than in the 
other cases with weaker fields, it does not scale with wave- 
length and remains almost at a constant value between 10 
and 100 ms~^ depending on filling factor. We interpret this 
behavior as a saturation effect in the sense that the line 
profile distortion due to a spot with B « 4000 G does not 
distort a measurement of the line center much more than 
a spot with B w 1000 G does. The reason for this is that 
Zeeman broadening already is so significant that the es- 
sential effect in the line profile is similar to a very cool 
spot in which the fiux from the spot area is simply miss- 
ing. For a field strength oi B — 4000 G, the displacement of 
a-components is wzocman 6kms^^ (Eq-H)), which means 
that flux from the spot area is essentially removed to the 
wings of the spectral line (note that the typical line width 
in Fig.[T]is ~ lOkms"^). In other words, the radial velocity 
signal does not grow any further as soon as the amplitude 
of Zeeman splitting from the spot area is comparable to the 
line width of the rotating star. 

The conclusion from this exercise is that the Zeeman sig- 
nal introduced by a magnetic spot that is not cooler than 
the rest of the star can be significant for radial velocity sur- 
veys aiming at precisions on the order of a meter per sec- 
ond. The RV signature has the same sign as the signature 
of a cool (non-magnetic) spot. The Zeeman signal grows 
with wavelength for moderate values of B inside the spot. 
The effect of a (non-magnetic) cool spot also can be very 
significant on the ms~^ level, but it scales with opposite 
sign, i.e., it is large at short wavelengths but diminishes to- 
wards longer. For a similar simulatio n using a toy model o f 
(non-magnetic) cool spots, we refer to lReiners et al.l (|2010t) . 
In their Fig. 10, they show that in a cool star (T ~ 4000- 
6000 K), the effect of a / = 2% spot that is AT = 200 K 
cooler than the surrounding is on the order of 10ms~^ at 
A — 550 nm, and the amplitude is a factor of 3 smaller at 
A — 1800 nm. Furthermore, they show that the RV am- 
plitude is larger if temperature contrast is larger (cooler 
spots in that example) , but dependence on wavelength also 
becomes a lot weak er for larger temperature contrast (see 
iReiners et al.ll2010L for more details). 

4. Symmetric line broadening 

In addition to the signature of co-rotating spots, the effect 
of variable average magnetic field distributed over the entire 
star can also be very interesting. This could, for example, 
approximate the effect of a magnetic cycle in an active star, 
or it can be caused by stellar rotation since the average field 
visible at a given epoch can differ from the one seen in other 
snapshots while the field distribution is rather symmetric. 
However, in our model description of the atomic and molec- 
ular Zeeman effect the pattern of Zeeman splitting is always 
symmetric (see Section[5]). Thus, for a geometrically sym- 
metric field distribution and single lines, no radial velocity 
shift can be introduced because the spectral line's barycen- 
ter always remains constant. On the other hand, the sym- 
metry of the appearance of several lines is broken as soon 
as line blending occurs. This can be seen in Fig. [31 Here, 
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Fig. 3. Example of apparent velocity shift in magnetically 
sensitive spectral lines. Black solid line: spectrum contain- 
ing four lines of molecular FeH without magnetic fields in 
a, T — 3700 K atmosphere; red dashed line: same spectral 
lines influenced by a magnetic field oi B = 1000 G. Inset 
shows auto-correlation function of the field-free case (black 
solid line) and cross-correlation between field-free spectrum 
with magnetic case spectrum (red dashed line). 

the difference between the two spectra is due to a change 
in the average field not restricted to a starspot. Due to 
line blending, the two spectra are now very different, and 
the barycenter of the merit-function (the cross-correlation 
function or the goodness-of-fit estimate) will be distorted 
introducing a spurious line shifiQ. Deviations from the cen- 
ter of the non-distorted line may be significant on the m s~^ 
level, which is much less than 1/1000 of the line width. 

Sign and amplitude of the RV shift depends on the sen- 
sitivity of the lines to the Zeeman effect but also on the 
amount of blending between different lines. Therefore, the 
apparent RV shift likely scatters between different wave- 
length bands and the net result is probably enhanced (ran- 
dom) jitter in the RV signal instead of a systematic RV 
offset like in the case of co-rotating spots. We calculate this 
effect and limitations predicting its amplitude in Sect. 16.31 
To do so, we require accurate spectral line data and a re- 
alistic description of the Zeeman effect over large ranges of 
the stellar spectrum. We describe the line data we use for 
our more sophisticated model in the next Section. 

5. Model atmospheres and line data 

For the detailed calculation of stellar spectral lines over 
large wavele ngth regions, we used t he MARC^ model at- 
mospheres ()Gustafsson et al.l I2008D using three different 
models: T^s = 5750 K, log(g) = 4.5 (a solar-type star), 
Teff = 3700 K, log(5) = 5.0 (an early-M-type star of 



The amplitude of this shift will in fact depend on the method 
the barycenter is located: fitting the dashed line in the inset 
in Fig. [3] will provide a different result than searching for the 
absolute maximum of that function. 
^ http://marcs.astro.uu.se 
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Fig. 4. Average Lande g factors for lines deeper than 2 % of the continuum in the modeled spectra range (blue circles) 
together with rms scatter of all Lande g values considered in that region (error bars). The three panels show three 
different atmospheric temperatures as indicated in the figure. The low values at 2.3 /xm are due to the CO band. 



spectral type -Ml), and Tcff = 2800 K, log(g) = 5.0 
(— M6.5 star). F or all models, solar abundances according to 
lAsplund et all (jlOQS) were assumed. Our theo retical spec- 
tra w ere computed using the SYNMAST code (jKochukhovl 
[200I that can treat atomic and molecular transitions in a 
magnetized medium. 

Atomic line data were ex t racted from the VALD 
database (jPiskunov et al.l 119951: iKupka et all Il999t) . The 
Lande g-factors were taken as provided by VALD or were 
computed using available term designations assuming the 
LS-couphng approximation. 

The spectra of cool stars are dominated by molec- 
ular absorption. However, for most molecular lines the 
Zeeman effect is only poorly understood. The lack of lab- 
oratory measurements of Lande g-factors and the complex 
physics behind molecular line formation in plasmas with 
strong magnetic fields make it very challenging to accu- 
rately model spectra of these objects. For this work, we con- 
centrate mostly on the effect of molecular FeH. We tested 
theoretically computed Zeeman patterns on obse rvations 
of a n umber of M-dwarfs as descri bed bv |Shulvak et al.l 
(|2010f) . FeH line data was taken from lDulick et al.l (|2003H 
an d we used corrected Hne intensities and positions follow- 
ing |^fende||eFal] (|2010l ). We also included FeH lines in the 
range 1.0-1.7/im usin g the same procedur e for computing 
Lande g-factors as in Shulvak et al.l ()2010[ ). In addition to 
FeH lines, we also included t he line list of CO transi- 

tions from lGoorvitchI (|1994[ ). The well-known 2.3/im band 
of CO is often used for RV m easurements beca use of its 
magnetic insensitivity (see, e.g. iBean et al.ll201flh . 

We did not consider Zeeman splitting of molecular 
bands other than from FeH and CO. Thus, the present in- 
vestigation is only an approximation of the effect of Zeeman 
splitting on RV signals; a complete model should include all 
molecular lines present in the stellar spectra. Nevertheless, 
FeH is the most important opacity source in near-infrared 
spectra of very cool stars, and we believe that the main 
effects from Zeeman splitting can be captured by our ap- 
proach. Prominent molecular bands not included in our 
model are due to, e.g., TiO, CH, OH, and MgH. Some of 
them are known to exhibit moderate or strong magneti c 
sensitivity as discussed in iBerdvueina fc Solankil (j2002l) . 
Since our main goal is to characterize the general behavior 
of RV signals comparing different wavelength regions, and 



http:/ /bernath.uwaterloo.ca/FeH/ 



since a large part of the trends can be described neglecting 
detailed line list information (see Sect.]!]), we do not see a 
reason why the addition of more molecular species should 
systematically change our results. Nevertheless, quantita- 
tive predictions about absolute RV distortions need to be 
interpreted with great care bearing in mind the limits of our 
modeling approach. Including more magnetically sensitive 
lines will also result in stronger blending, hence implications 
from line blending presented in this work are probably lower 
limits. 

Fig. U shows a compilation of average Lande g values 
used in our model spectra. For parts of 100 nm length, we 
calculated average Lande g values taking into account all 
atomic and molecular lines stronger than 2 % (neglecting 
line broadening due to surface motion and stellar rotation) . 
The figure also visualizes the spectral regions we considered 
for our investigation, which are similar to the photometric 
bands V, /, Y, J, 7J, and K. Cooler stars show more molecu- 
lar absorption than hotter stars do. This becomes apparent 
in the lower values of average Lande g values since on av- 
erage molecular lines have lower g. 

6. Results 

The mechanism causing an RV signal in an idealized spec- 
tral line through the Zeeman effect was demonstrated in 
Sect. [3] The amplitude of an RV signal occurring in a real- 
istic spectrum will also depend on blending and the prop- 
erties of lines contributing to radial velocity information, 
foremost their Lande factors and central wavelength. The 
influence on radial velocity measurements in real stars can 
therefore be expected to differ significantly from the results 
of our simple toy model. Furthermore, radial velocity shifts 
from temperature contrast and from Zeeman splitting lead 
to amplitudes that are comparable in absolute values if cal- 
culated independently, but active regions on the Sun are 
magnetic and differ in temperature. In analogy to the solar 
case, we expect that the stronger magnetic fields on other 
stars are present in spots that are both cool and magnetic, 
but it is very difficult to assess realistic values for tempera- 
ture and magnetic field contrasts. In the following, we first 
calculate the radial velocity signature of a magnetic spot us- 
ing not a single line but synthetic spectra including several 
thousand atomic and molecular lines that are split due to 
the Zeeman effect. After that, we show a few example cases 
for magnetic, cool spots and the net radial velocity signal 
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Fig. 5. Zeeman RV signature of co- rotating magnetic spots 
calculated from polarized radiative transfer over several 
wavelength bands. Stars and filled circles indicate six wave- 
length bands where the RV signature was measured, a fit 
to the first five bands using a straight line is shown as a 
solid line for each case (the CO band is not included in the 
fit). Dashed lines show the scaling expected from our our 
toy model, i.e., a = 2 in Eq.|3l 



from the two competing effects, and we calculate the RV 
signal of average magnetic field variations using synthetic 
model spectra. 

6.1. Co-rotating magnetic spot 

We calculated model spectra of rotating stars with mag- 
netic spots as in our toy model above, but for the more 
realistic case we used spectra from our polarized radiative 
transfer code. As above, the cross-correlation function is 
calculated using the non-spotted (non-magnetic) star as a 
template and the star with a magnetic spot as our data 
set. For each case, we calculated the bary center of the 
cross-correlation function in six wavelength areas of more 
than 100 nm each (see Fig.|3]). Three stars of temperature 
T = 5750 K, T = 3700 K, and T = 2800 K are calcu- 



lated, spot sizes of / = 1 %, 5 %, 10 %, and 20 %, and field 
strengths inside the spot oi B — 100 G and 1000 G are used. 
The results are shown in Fig.[S] 

For our hottest example, the results from the synthetic 
model do not differ very much from the trends we found 
in our toy model using a single line with constant Lande 
g. In all cases, spots with magnetic fields of 1 kG introduce 
RV signals exceeding lms~^ at all wavelengths. The scal- 
ing of Eq.[5] is a very good approximation to the situation 
seen in the case for T = 5700 K. In our cooler examples, 
however, the values we calculate from the synthetic atmo- 
sphere models show significantly less dependence on A in 
the sense that the radial velocity signal does not grow pro- 
portional to A squared. Nevertheless, we can approximate 
the radial velocity signal in all three examples using the 
following formula: 

A«..d = 300^/ (^X(—X, (3) 
s V kG / \ /im J 

with the relative fraction of the spot area, /, the magnetic 
field inside the spot, B, and the wavelength A. The value of 
the constant results from the geometry of the surface and 
on the distribution of Lande factors across the wavelength 
range, but not on other free parameters. While a « 2 in 
our hottest example (same value as in our toy model), we 
find a < 2 for cooler stars. The reason for this is clear 
from the distribution of Lande g-values shown in Fig.|4l 
While in the sun-like case (T = 5750 K) the typical Lande 
g-values are not a function of wavelength (the CO band 
being the only exception), Lande g-values are significantly 
lower at longer wavelengths in the cooler stars, partially 
compensating for the linear increase of the total Zeeman 
shift as a function of wavelength (see Eq.[T|). This is because 
at cooler temperatures, molecular lines become more and 
more important in relation to atomic lines, and (at least in 
our model) the molecular lines tend to have lower Lande 
g-values on average. 

A central result of our simulation is that the radial ve- 
locity signal due to a co-rotating spot of 1 % the size of the 
projected stellar disk and a magnetic field strength of 1 kG 
has an amplitude of approximately 3ms^^ if observed at 
A = 1 /im. The amplitude scales linearly with filling factor 
and quadratically with both magnetic field strength and 
wavelength. In stars significantly cooler than the Sun, the 
scaling with wavelength is weaker than quadratic because 
the relevant spectral features are magnetically less sensi- 
tive. 

6.2. Cool and magnetic spot 

In the preceeding sections and in iReiners et al.l (|201flh . 
the effects of magnetic spots (at the same temperature as 
the photosphere) and of (non-magnetic) temperature spots 
were modeled independently. In reality, of course, spots are 
understood to be both magnetic and effectively cooler than 
the quiet photosphere. As a first approximation, we com- 
pared both effects. We found that temperature and mag- 
netic effects can cause radial velocity signatures of compa- 
rable amplitude, but temperature effects are larger at short 
wavelengths while magnetic infiuence is more significant at 
long wavelengths. 

The first order approximation looking at both effects in- 
dependently is probably not realistic because a low temper- 
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Fig. 6. Radial velocity amplitude due to co-rotating cool, 
magnetic spots (/ — 1%, B — 1000 G, v sini — 
2kms^^) using our synthetic atmosphere model for differ- 
ent star/spot temperatures. Black circles show the influence 
of temperature only (_B = 0), red circles show the net effect 
including temperature contrast and Zeeman broadening in 
the spot. Blue crosses are average values and range for six 
wavelength bands. Green triangles show residuals between 
red and black circles, i.e., the effect due to Zeeman broaden- 
ing. Left panel shows a sun-like star with T = 5750 K, right 
panel shows an early M-type star with T = 3700 K. Top to 
bottom panels show different values of AT = T^tar — ^spot; 
top panel: AT = 200 K; middle panel: AT « T/5; bottom 
panel: AT « T/3. 

ature spot emits less flux than a hotter one, which will lead 
to a weaker influence of the Zeeman effect on the line pro- 
file and its dependence on wavelength. In order to consider 
the two effects in a consistent way, we calculated models 
of stars with spots that are both cool and magnetic, and 
we show the results of radial velocity amplitude in Fig.|6l 
In all cases, spot parameters are f — 1% and B = 1000 G, 
the star is assumed to be rotating at v sini = 2kms~^. We 
show results for a sun-like star with T = 5750 K and for a 
cooler (M-type) star with T = 3700 K; the T = 2800 K does 
not provide new information and we do not include it in this 
example. For spot temperatures, three cases are considered 
for each star: one with AT = Tgtar — Tspot = 200 K (top 
panel), one with AT w T/5 (middle panel), and a third 
with AT w T/3 (bottom panel). 

Our first result is that in all cases the net effect is a 
non-trivial combination of temperature contrast and mag- 
netic Zeeman splitting. Both mechanisms work in the same 
direction; the total RV signal accounting for both effects is 



always larger than the signal from one mechanism alone. 
The net effect reaches values up to the 10ms~^ level. In 
our examples with the lowest spot/star temperature ra- 
tio (AT = 200 K), the radial velocity signal monotonically 
grows with wavelength. In our intermediate cases of tem- 
perature contrast (AT « T/5, center row), the net sig- 
nal has a local minimum around 1 /im, it is dominated by 
temperature contrast effects at shorter wavelengths and by 
Zeeman splitting at longer wavelengths. In our examples 
with the highest temperature contrast (bottom panel), the 
wavelength dependence of the total RV signal is nearly con- 
stant with wavelength: in the sun- like star (bottom left), 
the RV signals from the temperature effect and from the 
Zeeman effect almost cancel; in the cooler star with the 
coolest spot (bottom right), the RV signal is always domi- 
nated by the temperature contrast. In the latter case, the 
spot contributes so little flux that the Zeeman signal is 
not significant even at the longest wavelengths used here. 
Interestingly, the signal from the temperature spot alone 
shows some scatter with wavelength that is not monotonic 
in A and is produced by the temperature dependence of 
individual absorption lines that sho w different intensitie s 
inside and outside the spot region (cp. iReiners et alll2010[) . 
For example, the CO lines in this case become deeper with 
lower temperature, which counteracts missing line absorp- 
tion emerging from the spot area and leads to a reduced 
RV signal within the CO line region. 

Even if the model we show in Fig.[S] includes radiative 
transfer of atomic and molecular spectral lines, and both 
the inffuence of temperature and Zeeman splitting, we are 
aware of the problem that our model is probably very dif- 
ferent from any real star. Values of spot temperature and 
field strengths in our model are probably not unrealistic 
for some individual spots, but the real range in tempera- 
ture and field strengths are largely unkown. Perhaps more 
important, active stars most probably are very different to 
a one-spot model and evolve in time. Therefore, we restrict 
ourselves to the examples shown in this section and do not 
attempt to make any more specific predictions on radial 
velocity signals in active stars. 

We conclude from our examples that in active stars, 
both the influence of temperature contrast and the Zeeman 
effect can be of comparable amplitude, and that the de- 
tails of spot distribution and their temperatures and mag- 
netic fields determine amplitude and wavelength-scaling of 
the RV signature due to activity. Thus, radial velocity sig- 
nals due to active regions cannot be expected to vanish at 
infrared wavelengths. On the contrary, in many stars the 
infiuence of starspots may be much more severe at longer 
wavelengths. The best way to discriminate between a planet 
and a starspot is therefore simultaneous measurement of ra- 
dial velocities at many different wavelengths; even if the 
scaling of the radial velocities with wavelengths is diffi- 
cult to predict, it is improbable that a signal due to co- 
rotating active regions is independent of wavelength. Any 
wavelength-dependent signal will rule out companions as 
source, and the scaling with wavelengths will provide use- 
ful information on the nature of active regions. 

6.3. Average field variations 

Magnetic and cool spots on the surface of a rotating star 
introduce line profile variations as discussed above. Radial 
velocity signals due to this effect are introduced by the de- 
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Fig. 7. Radial velocity scatter due to variable average magnetic field. Red circles show results from cross-correlation 
between spectra of a non-magnetic star and a star with, an average magnetic field. Blue crosses and error bars show 
average and rms-scatter of individual chunks for six wavelength bands. Left panel: surface fraction / =10% of the star 
is covered with magnetic field varying between and 100 G; right panel: the entire star (/ =100%) shows magnetic field 
variability of 100 G. The small scatter at 2.3 fim is due to the CO-band. Note the different scaling on the y-axis. 
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Fig. 8. Zeeman-induced RV scatter for the case of variable 
average magnetic fields. Six cases are shown in which the 
fractional coverage of the field is between 1 % and 100 %. In 
that area, the field varies by 100 G. The scatter depends on 
wavelength and can be approximated by the dashed grey 
lines calculated from Eq.m 



formation of all individual spectral lines and show the same 
period as the rotation of the star. Another effect caused by 
the Zeeman effect was introduced in Section|4] and is due 
to systematic mismatch between large spectral regions ob- 
served at one time with respect to another observation. A 
possible reason for such a mismatch can be variation in 
the (average) magnetic field of a star, for example during 
a magnetic cycle. In a single line, average-field variability 
would not lead to a radial velocity signal assuming symmet- 
ric Zeeman splitting (as always assumed here). In spectral 
regions containing many lines, however, an apparent shift 
may be introduced because lines are usually blended with 
others (see Fig. [3]). 

In Fig. [71 we show the results from cross-correlating a 
spectrum with no magnetic field with a spectrum of aver- 
age field of Bf = 10 G (left panel) and Bf = 100 G (right 
panel). For field variations with respect to a non-zero field 
we expect similar results since Zeeman splitting is linear in 
B. The spectra are constructed assuming B = 100 G field 
strength in active regions homogeneously covering / = 10 % 
(10 G case) and / =100 % (100 G case) of the projected stel- 
lar disk. The examples represent the cases in which the aver- 



age magnetic field varies by AB = lOG and AB = 100 G, 
respectively. We show results from individual wavelength 
regions covering 10 nm each. In contrast to the cases with 
co-rotating spots, the radial velocity shift does not follow 
a systematic pattern because the signal we measure here is 
a result of the random blending of broadened lines. On av- 
erage, the signal is zero but variations in the field strength 
introduce scatter that depends on the properties of the field 
variations and wavelength. In Fig. [3 we show average val- 
ues (blue crosses) and rms-scatter for wavelength bands 
containing several adjacent individual wavelength parts. As 
expected, we find that the average values of the RV shift is 
consistent with zero at all bands, but the rms grows with 
wavelength (an exception is again the CO band at 2.3 /zm). 

In Fig. [51 the rms- values of the RV signal are displayed 
as a function of wavelength for the two models of Fig.[71plus 
four more models with fractional coverage between / =1 % 
and 100%. Note that the amplitude of rms scatter is given 
in cms~^. We fit a power law to our results (grey dashed 
lines in Fig. [8]) and find an approximation to the RV scatter 
introduced by variations in the average magnetic field: 

/ \ \ 118 

CT(«rad)-10— / — , (4) 

s \ /im / 

for AB = 100 G. This approximation is valid for all wave- 
length bands except the one containing the CO lines. For 
the AB = lOG case, we find that while the effect is of 
course smaller than the one for AB = 100 G, there is no 
simple scaling relation that relates the cases of different 
AB. Nevertheless, we can conclude that additional scat- 
ter in radial velocity measurements can be introduced by 
variable average magnetic fields. The uncertainty of a radial 
velocity measurement at a given wavelength can be affected 
by field variability, but the amplitude of this effect is likely 
well below lms~^ in realistic cases. Even stars with very 
strong magnetic fields are not believed to show variations in 
average field strength on the order of kilo-Gauss, and such 
stars would probably show very strong variability in chro- 
mospheric emission lines, too, which could help to identify 
such cases. 

7. An active star example: AD Leo 

Our calculations predict that magnetically active stars may 
show wavelength-dependent RV variations with a larger RV 
signal at longer wavelengths. For a first test of our results. 
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we searched the HARPS-ESO arehive for magnetically ac- 
tive stars with clear periodicities detected on them. A very 
clear and prominent example we found was AD Leo (Gl 388, 
Spectral Type M4.5Ve). HARPS observations on this star 
span a time baseline of 900 d and have a typica l SNR of 
50 at 600 nm. As reported bv lBonfils erall (|201lD . AD Leo 
shows a very strong periodic signal in the RVs at 2.22 d, 
which is consistent with the stellar rotation period found 
in the Zeeman Doppler Imaging analysis by iMorin et al.l 
(|2008D . At the present time, 40 spectra are publicly avail- 
able. We analyzed the data using t he HA RPS-TERRA soft- 
ware (jAnglada- Escude fc Butler! |2012[ ). HARPS-TERRA 
derives RV measurements by least-squares matching each 
observation to a high signal-to-noise (SNR) template gen- 
erated by CO- adding all available observations. It also allows 
to obtain the Doppler measurement using a limited num- 
ber of echelle orders at a time, enabling the analysis of RV 
signals as a function of wavelength. In order to accumulate 
enough SNR to derive good quality RV measurements, we 
split the stellar spectrum in 7 parts using 419, 450, 486, 
528, 583, 631 and 665 nm as the central wavelength of each 
part. Except for the redder two parts, each part spans ten 
HARPS Echelle orders, that is: orders 10-19, 20-29, 30-39, 
50-59, 60-66, 68-71. The last two parts are chosen to avoid 
order 67 containing Ha that is highly variable in active 
stars. In all these par ts, the 2.22 d pe riod is clearly detected 
in the periodograms (jCummin 3 l2004[ ). The 10 bluest echelle 
apertures (0-9) are not discussed here because the SNR is 
very low at the bluest orders (at the 5th echelle aperture 
it was typically below 5), and uncertainties associated to 
each individual RV were of the order of 50ms~^. 

In the upper panel of Fig.|9l we show our fit with period 
P=2. 22704 d to the RVs derived for each part. For simplic- 
ity, a sinusoid with the fixed period derived from the best fit 
to the RVs using the full spectrum (see top periodogram in 
Fig.[Tni) was adjusted to each part, so the only free param- 
eters were amplitude and phase of the signal. In the lower 
panel of Fig. [HI we plot the derived semi- amplitudes as func- 
tion of central wavelengths for each part. The uncertainties 
were derived using the bootstrap technique, i.e., comput- 
ing the scatter of the amplitude as obtained by randomly 
selecting samples of the observations. The HARPS-Data 
Reduction Software also provides a measure of the mean 
spectral line asymmetry, called the bisector span (BIS). 
BIS is a measure of the asymmetry of the cross-correlation 
function in RV space as obtained from cross-correlating the 
stellar spectrum to a binary mask (M2 binary mask, see 
Pepe et al. 20021. As demonstrated in several studies (e.g. 
Queloz et al. 2001) , BIS often anti-correlates with spurious 
RV offsets if a cool spot is responsible for the apparent 
RV shifts. The BIS periodogram of AD Leo shows 4 peaks 
of similar power at 1.813, 1.950, 2.041, and 2.219 d (bot- 
tom panel in Fig. lTUl) . These periods are all related through 
daily aliases (standard and sidereal day) and they likely 
correspond to the same physical periodicity. Since none of 
the BIS periods exhibits a false alarm probability lower 
than 1%, the periods are subject to significant uncertain- 
ties; only an approximate match can be done when com- 
paring the BIS measurements with those in the RV signals. 

The photometric period (and, presumably, the rota- 
tion period of AD Leo) has been reported to be 2.23 d 
(|Engle et al.l I2009D favoring the period of 2.219 d as the 
most likely fundamental signal in the BIS. To see how it 
compares to the RV signal, we fixed the period of a sinu- 
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Fig. 9. Upper panel: Phase folded radial velocities and fit- 
ted signals to the preferred Doppler period of 2.22704 d. 
The phase folded fit to the BIS is also provided. For illus- 
tration purposes the BIS data has been shifted and multi- 
plied by 10. Right panel: Semi-amplitude K of the signal 
as a function of wavelength. The best fit to a linear model 
representing the wavelength dependence of K is given as 
a thick dashed line. The semi- amplitudes derived from the 
full spectrum (using CCF and HARPS-TERRA measure- 
ments) are also show as solid brown circles. Note that both 
measurements provide incompatible results providing a fur- 
ther test to assess the reality of a Doppler signal. 



soid to the preferred RV period and adjusted the ampli- 
tude and phase of the BIS curve (see Fig. [51). Doing this we 
found that BIS appears to be anti-correlated with the RV 
curve, which is consistent with the expectations for a cool 
spot-induced signal. Note that, for visualization purposes, 
the BIS measurements in Fig. [9] were multiplied by a fac- 
tor of 10 (and are offset). Also note that while the 2.22 d 
periodicity is clearly detected in the RVs, the F-ratios of 
the BIS candidate signals are 15 times smaller and barely 
significant (see Fig.[Tni). We also looked at other activity 
indicators typically associated to spurious RV signals on 
M dwarfs (e.g., FWHM of the CCF or the S-index; see 
iLovis et al.ll201lHAnglada-Escude et al.ll2012[ ). but we did 
not find any further indication of a peak near 2.22 d. 
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Fig. 10. Top: Periodogram of the radial velocities ob- 
tained using the full spectral range of HARPS. Bottom: 
Periodogram of the BIS. The 1% FAP threshold is illus- 
trated as solid horizontal lines in both panels. The four 
peaks in the bottom panel are likely strong aliases of the 
same signal (2.22 days would be compatible with the pho- 
tometric period of the star). Note that their average FAP 
is only 5%. 



Assuming that the Doppler signal is caused by cool 
spots (neglecting Zeeman broadening), one would expect 
its amplitude to become weaker at redder wavelengths. We 
obtained a weighted fit to our values of K shown in the 
bottom panel of Fig. [3] using a linear model of the form 
K[X] = AX + b, where A is the slope and is measured in 
units of velocity per unit of wavelength (ms~^ /Ltm~^). The 
obtained slope is 26.4 m s~^ /i m^^ which is positive and dif- 
ferent from at a level. Therefore, we conclude that the 
RV signal does not diminish towards longer wavelengths in 
the range covered by HARPS. Instead, the obtained wave- 
length dependence of the amplitude has a significantly pos- 
itive slope with larger amplitudes K towards longer wave- 
lengths. Our data do not extend into the infrared wave- 
length range, and we cannot reach any firm conclusion on 
the behaviour of the RV signal wavelengths longer than 
700 nm. Nevertheless, our example provides first evidence 
that the RV signal of an active star does not always dimin- 
ish at larger wavelengths. Since RV signals from cool spots 
due to temperature contrast alone are supposed to show a 
monotonic behavior across optical and infrared wavelength 
regions, we interpret this as evidence for another mecha- 
nism causing the positive slope of the RV curve. The results 
of this paper would indicate that the Zeeman effect is the 
reason for this, and that the amplitude of this signal would 
be even larger at infrared wavelengths. 

We note that the semi-amplitude derived from CCF 
measurements differs from the one derived from HARPS- 
TERRA RVs using the full spectrum; the result from CCF 
is K ^ 30.5 ± 1.0ms-\ the resuU from HARPS-TERRA 



is K = 23.4 ± 0.9ms^^. This indicates that the changes 
in the line shapes affect each method in a very different 
way, further indicating that the measured RV offsets are 
due to changes in the line profile shapes rather than real 
Keplerian signals. In conclusion, even in this example, in 
which we know the photometric period and see an anti- 
correlation between the RV with line asymmetries, the as- 
sumption that the RV signal is induced by the tempera- 
ture contrast effects alone is likely to produce an incorrect 
interpretation of the data and can lead to the erroneous 
prediction that the spurious RV signal is suppressed at nIR 
wavelengths. This misconception would have serious conse- 
quences if one attempts to correct the RV curve for activity 
signals. Nevertheless, it is certainly true that Keplerian sig- 
nals cannot be wavelength-dependent. If stellar activity on 
a time-scale similar of a Doppler signal is suspected, only 
a comprehensive analysis of its wavelength dependence can 
shed light on its true physical origin. 

8. Discussion 

We have investigated the influence of magnetic activity on 
radial velocity measurements in active stars. In contrast 
to earlier calculations, we included the Zeeman effect and 
calculated line barycenter shifts due to spots that are cool 
and magnetic. Furthermore, we looked at the case of vary- 
ing average magnetic fields that are not concentrated in 
co-rotating regions and may introduce signals not in phase 
with stellar rotation. 

Our most important result is that co-rotating magnetic 
starspots can be expected to significantly distort stellar line 
profiles and RV measurements. Neglecting the cool temper- 
ature of the spots, the signal from the Zeeman effect alone 
easily exceeds the lms~^ level even for very small active 
regions (/ = 1%) and slow rotation {v sini = 2kms~^) 
if the field inside the spot is comparable to sunspot fields 
{B ^ IkG). This signal has the same sign as the signal 
from a non-magnetic, cool spot, it grows with both mag- 
netic field and wavelength, and is approximately four times 
larger at A = 1000 nm than at A = 500 nm. Its amplitude 
saturates for very strong fields above approximately 1 kG. 
Comparable radial velocity signals are also expected from 
the temperature contrast of cool spots alone (neglecting the 
Zeeman effect) as shown in earlier investigations, but they 
are largest at short wavelengths. We note that the effect of 
Zeeman broadening on integrated (non-polarized) light is 
independent of magnetic polarity. Therefore, the RV signa- 
ture does not depend on magnetic polarity, and in particu- 
lar does not cancel out if magnetic areas consist of several 
spots with opposite polarity. 

In contrast to the systematic fake RV signals from co- 
rotating spots, the RV signal due to variable average (non- 
localized) fields is of statistical nature only affecting the 
noise floor of the measurements. The typical uncertainty of 
jitter induced by an average magnetic field varying by 100 G 
is on the lOcms""'^ level. No systematic signal is expected 
from variations in the mean magnetic field strength. 

Including both temperature contrast and the Zeeman 
effect is necessary to understand RV signals in active stars. 
In case of large contrast (very cool spots), the Zeeman effect 
is less important because the integrated spectrum contains 
less fiux from the magnetic (spot) region. Nevertheless, 
since temperature contrast diminishes at long wavelengths 
and the infiuence of Zeeman broadening grows in the same 
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direction, the RV signature of the Zeeman effect is signifi- 
cant in infrared observations of active stars. Thus, the RV 
signal of active stars does not vanish at long wavelengths, 
and infrared observations are not less affected by activity 
than observations at optical wavelengths. This may be par- 
ticularly important for moderately active stars that could 
be p opulated by magnetic areas in a way similar to the Sun 
(jLagrange et al. 201Q,) . 

The magnetically insensitive CO lines in the K-band 
provide a notable exception. Their response to magnetic 
fields is so low that even strong fields do not substantially 
distort their line profiles. These CO lines can be very useful 
to disentangle RV variations due to Keplerian orbits from 
magnetic activity. Unfortunately, the CO lines are contam- 
inated by telluric lines introducing oth er problems when an 
accuracy on the ms~^ level is desired (jBean et al.ll201Cll ). 

Radial velocity signatures due to conve ctive blueshift 
can also be significant in sun-like stars ()Meunier et al.l 
l201Clt) . Their amplitude likely depends on line-depth and 
therefore adds additional complexity to disentangling stel- 
lar activity from Keplerian signals. 

It is very difhcult to predict the dependence of the RV 
signal as a function of wavelength because it sensitively de- 
pends on the combination of spot temperatures and their 
magnetic field strengths. Both are poorly constrained by 
currently available data and simulations suggest that differ- 
ences between sola.r and very cool star magnetic structures 
exist (|Beeck et al.l l201l[ ) . Simultaneous measurements of 
RV amplitude over large wavelength regions provides useful 
information for characterizing stellar activity, most impor- 
tant starspot temperature and magnetic fields. In moder- 
ately active stars, the precision required for such a measure- 
ment is on the level of a few ms~^ for wavelength intervals 
of several hundred nm, which is a challenge for typical RV 
surveys. Simultaneous RV measurements at different wave- 
length bands are possible already in a few spectrographs 
and will become accessible over very large ranges with high- 
precision RV spectrographs operating at infrared wave- 
lengt hs, as for exa niple CAR MENES (jOuirrenbach et al] 
l2010t) and SPIRou (jArtigau et al.ll20lTI) . Data from these 
facilities will not only provide a reliable method to distin- 
guish a Keplerian signal from magnetic activity, they will 
also allow a deep look into the magnetic and temperature 
structure of stellar surfaces. 
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